Non-equilibrium interface equations: An application to 
thermo-capillary motion in binary systems 



Ravi Bhagavatula^ , David Jasnow^ and T. Ohta* 
''' Department of Physics and Astronomy, University of Pittsburgh 

Pittsburgh, PA 15260 
Department of Physics, Ochanomizu University, Tokyo JAPAN 112 



February 1, 2008 



Interface equations are derived for both binary diffusive and binary fluid systems subjected to non- 
equilibrium conditions, starting from the coarse-grained (mesoscopic) models. The equations are 
used to describe thermo-capillary motion of a droplet in both purely diffusive and fluid cases, and 
the results are compared with numerical simulations. A mesoscopic chemical potential shift, owing 
to the temperature gradient, and associated mesoscopic corrections involved in droplet motion are 
elucidated. 



keywords: interface equations, thermocapillary motion, capillarity 



1 



I. INTRODUCTION 



The study of interfacial dynamics in multi-phase systems has been a topic of considerable interest in recent years. 
A wide range of phenomena such as solidification, viscous fingering, droplet migration, spinodal decomposition and 
fracture all involve, in one way or another, a description of interfacial motion. Often analysis of these and related 
phenomena is performed at the macroscopic level, using, for example, hydrodynamic equations supplemented by 
appropriate phenomenological boundary conditions. Macroscopic interface equations are exceptionally useful both 
conceptually and computationally in that bulk degrees of freedom are eliminated; computationally the cost is the 
required careful tracking of the interface position and shape. Some examples of their application are included in 
Ref. [Q and citations therein. 

The macroscopic approach, while valuable in its own domain, has natural limitations. The most obvious lie in 
situations in which the scale of structures is not sufficiently larger than the physical interfacial thickness, which can 
grow to several molecular length scales not too far from criticality, for example. Phenomena such as coalescence 
ultimately run macroscopic approaches to their limits, although ad hoc but reasonable methods are often available for 
resolving singularities having to do with infinitely sharp interfaces. In addition, one may wonder whether macroscopic 
boundary conditions may require correction in situations which strain the assumptions. 

Derivations of macroscopic interface equations from more microscopic starting points can be extremely useful in 
illuminating just what assumptions are necessary and the potential limitations. Indeed considerable progress has 
been made in this regard beginning from a coarse-grained or mesoscopic level. [|| However, the interfacial response 
to non-equilibrium perturbations or inhomogeneities such as an imposed small thermal gradient, can pose interest- 
ing challenges. In this paper, we develop interface equations starting from mesoscopic (coarse-grained) models for 
both purely diffusive binary systems and binary systems with hydrodynamic interactions, subjected to small non- 
equilibrium perturbations. Using this approach, we investigate as an example, thermo-capillary driven motion of a 
droplet in both binary diffusive and binary fluid systems in an imposed temperature gradient. Note that motion of a 
droplet in the hydrodynamic case has been well explored by the use of macroscopic techniques , while diffusion 
driven thermo-capillary motion has not been received much attention in the literature Thus, to the best of our 
knowledge, our modeling approach provides a first detailed comparison for the droplet motion in these two distinctly 
different cases. Furthermore, in addition to reproducing the results consistent with macroscopic hydrodynamics, the 
coarse-grained modeling approach also provides a means to investigate mesoscopic corrections associated with finite 
interfacial thickness. This feature may prove useful in understanding interfacial dynamics at mesoscopic length scales. 

The remainder of this paper is laid out as follows. In the next Section we introduce the notation and briefiy 
review the coarse-grained models we will consider along with the relevant equilibrium features. Section [II deals 



with the introduction of the types of non-equilibrium considered here and with the effect on a single flat interface. 
Section IV contains the derivation of the macroscopic interface equations and mesoscopic corrections in both the fluid 
and diffusive cases. 



II. MODELS 



We begin with the d-dimensional coarse-grained (mesoscopic) binary fluid and binary diffusive models which are 
commonly referred to as the Model-H and Model-B in the literature . These traditional models have been extensively 
used to investigate critical fluctuations in binary systems. However, the usefulness of such models for regions well below 
the critical point is still being established, by, for example, demonstrating that in particular limits well established 
macroscopic results are reproduced. |^ This remains a focus in this paper. One advantage of this type of modeling 
is that the physics of capillarity is naturally built-in as opposed to being introduced through boundary conditions in 
macroscopic techniques. 



A. Hydrodynamic Model 

The coarse-grained model is defined by a conserved order parameter {e.g., concentration of one of the phases) 
whose evolution is given by 



dt 



+ V ■{(bv)^\7'^iir,t) (1) 



where fi — (§^) is the appropriate chemical potential given by the functional derivative of the Helmholtz free energy 
F. The fiuid velocity t/(r, t) satisfies a modified Navier-Stokes equation, 
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p— = Tj\/^v- VP + /iV0 (2) 

where P is the pressure, rj and p are viscosity and density respectively, which are assumed to be fixed and equal for 
both phases for simplicity. Equations (|l|) and (|^) along with the incompressibility condition, V ■ v — 0, completely 
specify the hydrodynamic model with the following three boundary conditions: (Bl) dfi/dn — 0; (B2) d4)/dn — and, 
for example, (B3) v = Q on the boundary. Here d jdn corresponds to the normal derivative at the system boundary. 
The condition Bl preserves the global conservation of order parameter 0, B2 is the natural boundary condition which 
demands smoothness of (j> at the edges, and B3 enforces "no-slip," with other possibilities easily included. Note that 
our interest here is to address small Reynolds numbers (Re) and, hence, the nonlinear convective terms have been 
dropped in the velocity equation @). 



B. Diffusive Model 



Model-B is in turn specified by evolution equation 

VV (3) 



d(t){r,t) _ 2, 



dt 

with boundary conditions Bl and B2 mentioned above. One may view this as the diffusion dominated limit of Eq. (Q) 
where the fluid velocity is neglected. This is a version of Cahn-Hilliard model for binary diffusive systems used to 
study phase separation 

The mean-field level equilibrium properties of the two coexisting phases are characterized by the local Helmholtz 
free energy density which, for systems of interest, we take as a polynomial function 



f = j2-r (4) 



with the coefficients q;„ depending smoothly on physical fields such as the temperature chosen such that the system 
is within the two phase regime. Although it is not essential we will generally assume a scaling equation of state (see 
below). The equilibrium properties depend on the local temperature T and other field variables that determine the 
an- For example, taking a2 = — t, where r > depending linearly on the temperature [r — (T^ — T)/Tc with Tc being 
the mean field critical temperature], and — 1 with all other — 0, we have the usual double- well Landau free 
energy, which will serve as our "standard model." The (Helmholtz) free energy functional for this case is given by 

„M ^ + nm^r ^ - 4 + *l^'r, (.) 

where we assume the mesoscopic length £ is typically of the order of several molecular spacings some distance below the 
critical point. This length characterizes spatial variations of the order parameter and consequently sets the interfacial 
thickness. The chemical potential fi is given by 

^(r, t) = -ev^4> + ^ = ~fV^4> - r<^ + 0^ (6) 

Note that our standard model corresponds to a symmetric situation with equilibrium phase values (j)eq = ±v^! 
coexistence curve is symmetric and, for example, response functions and other bulk properties are identical in the 
coexisting phases. When the system is in equilibrium, a flat interface with normal along x has the well-known profile 
(exact within this mean-field level in an infinite system). 



[x) = (j)eqta.'ah[ — j^], (7) 



jX — Xi)^ 

where xq is the location of the interface and ^ — tl\/T is the thermal correlation length (see, e.g., Ref. [p|jTc|] 
for details of this analysis). Note that the profile takes values ±(/)eq at x = ±cxd respectively. In a finite system 
this serves as a good approximation for the profile so long as the system size is much larger than the interfacial 
thickness ^. We will assume more generally scaling behavior for the order parameter profile, (\>{x) = (/)eg/i(2-^), 
with h{y) ±1 as ?/ ^ ±oo. Such behavior will follow from a scaling equation of state determined by /(0) pT[ , p^ , 
or quite generally with the identification of ^ with i. The profile given in Eq. (0) satisfies Eq. (^) with constant [i. 



3 



Note that /i = for a symmetric system at two-phase coexistence and for the above interfacial profile. Using these 
features, and a scahng form for the profile, it is straightforward to show |^,@ that the surface tension a is given 

by (7(r) = £^ J{^)^dx = fifj^^ where ei is a non-universal constant depending on details of the profile h{x) and 

taking the value 2\/2/3 in the standard 0'*-case. Note that the surface tension depends on the temperature (and 
ultimately on the local temperature) through the (peg and ^. Also, note that the analysis can be extended by choosing 
the coefHcients a„ such that the macroscopic limit (i.e., formally, £ — > 0) of ct exists. Physically the macroscopic limit 
obtains when the scale of any structures greatly exceeds the interfacial thickness, which is of order ^ ~ £. Specifically, 
for the above (/)^-symmetric model, the choice a2 = tq/£'^ and 04 = uq/I'^ (with uq = 1 for simplicity) ensures the 
finiteness of a as well as the bulk order parameter values in the macroscopic or sharp interface limit, ^ — s- 0. For 
specificity we will use the symmetric cj)'^ model described above. The first equality in Eq. ^) is general; we will use it 
and indicate below which features are model dependent. 

III. INTRODUCING DRIVING TERMS 

When the an are spatially uniform and chosen such that the system is on the bulk phase boundary, the system can 
evolve to a simple equilibrium state with a single homogeneous phase or a two-phase equilibrium with flat interface and 
volume fractions chosen according to the overall order parameter. Considerable progress has been made in the study 
of kinetics of the approach to equilibrium and of phase separation in such cases; see, e.g.,Refs. HJl^. Macroscopic 
interface equations have been derived from the coarse grained level in order to to explore the physics of phase ordering 
kinetics [p|Jl^]. However, when one or more of the a„ has a small spatial variation, the situation is somewhat different. 

In such a situation, the system may be able to reach an inhomogeneous equilibrium state depending on the nature 
of the variation and the boundary conditions. For example, in a closed system, a fluid in a gravitational field reaches 
an equilibrium state with a spatially varying density. However, in a sufficiently large system, such an inhomogeneity 
can yield interesting dynamical evolution long before the walls come into play, [ p^ An example which produces quasi- 
stationary states is that of capillarity driven droplet motion in a temperature gradient. ||8|]. Here, in the simplest 
case, a droplet drifts toward the 'hot' side with approximately steady state motion until reaching the vicinity of 
the system boundary. In this approximate steady state, the system is spatially non-uniform, but close to local 
equilibrium everywhere. In such cases, the interfacial dynamics are influenced by the spatial variation in the surface 
tension resulting from the spatial dependence of the q;„. Our interest here is to address this spatially inhomogeneous 
situation and present an approach for deriving reduced equations that describe the interfacial dynamics. 

For specificity we consider the following spatial dependence in the standard example mentioned above: t{x) — 
— Q!2(a^) = To — /3a; > 0, a4 = 1, and all other a„ = in the free energy given by Eq. (Q). This situation corresponds to 
imposing a temperature gradient on the system along x direction with the two ends kept at temperatures corresponding 
to To and tq — (3Lx where Lx is the size of the system along x. We restrict ourselves to temperatures slowly varying 
over the correlation length so that << 1. This is easily satisfied experimentally, and our analysis remains valid 
even for temperature gradients large by experimental standards. |^ Furthermore we consider only the situation in 
which the entire system is in the two-phase region. Thermocapillary phenomena associated with such type of thermal 
gradients have been numerically investigated recently in Refs. [p|, p^ , p"7P ] . Note that the spatial variation and hence, 
the temperature field, is fixed, corresponding to large thermal conductivity of the medium. Generalizations to treating 
the temperature field as an active dynamical variable are possible and will be treated elsewhere. 

Before considering the evolution of a complicated interface, it is useful to examine the properties of a stationary flat 
interface with its normal along x in a thermal gradient. In such steady state, which exists for the case of conserved 
order parameter, d(j)/dt = and v — 0. Hence following from Eq. (0), fi should be a constant since d^/dn = at the 
system boundaries. However, the steady state profile 0(x) (consistent with a non- vanishing, constant value of /i) is a 
combination of two ramps as shown from direct simulations of the coarse grained model in Fig. 1. Below we will use 
a local equilibrium approximant which describes the ramp profile very well. The location a;o of the interface depends 
on the initial order parameter content in the system. It is important to recognize that for conserved order parameter 
the profile of Fig. 1 represents an equilibrium state, while for nonconserved order parameter the 'kink' quite generally 
moves with constant velocity toward the hot side as order parameter is converted in such a way as to reduce the total 
free energy. |]l8| 

For the present (conserved) case, by assuming that the interface at a; = is at local equilibrium corresponding to 
the temperature t{x = 0), one can obtain fi to first order in the gradient /3 as follows. First, to simplify the algebra 
we rewrite t{x) — m(x)^ with m{x) — mo(l — bx). To 0{l3) there is no change if we identify b — f3/2 and mp ~ 1. 

For simplicity we set ^ = 1 and have — —j;^ — m'^{x)(/) + 0'^. Following Ref. (isj we remove the local equilibrium 
order parameter and make a non-linear transformation defining (t>{x) = m{x)r](z{x)). We set the function z{x) by 
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requiring that the resulting differential equation for r]{z) has unity for the coefficient of cPrj/dz^ . One finds, then, 
that z{x) = mo (a; — 6x^/2) and 



-^ = -'^-^l + V^ + {ihmo/m^)^ (8) 

This procedure effectively isolates the explicit dependence on the temperature gradient b (x (3. Now one can perform 
an order-by-order analysis, writing rj = ??o(^) + br]i(z) + . . . and = fio + bfii -|- . . .. At 0{b°) we have 

with, as expected, 770 — tanh(0/-\/2). The function 0o = TO(a;)77o(-z) corresponds to the local equilibrium ramp and is 
an excellent numerical approximation. In fact, if plotted on the same scale as the simulation result in Fig. 1, there 
would be virtually no visible difference. 



At 0(5) we have 



= Cr]i + — (10) 

TOn Too dz 



where the operator C — —d'^/dz^ — 1 + Stjq is recognized as the fluctuation operator for the 0"* theory, and drj^/dz is 
the so-called 'translation mode' (see, e.g., Ref. |l^). Multiplying through by drjo/dz, and recalling that Cdrjo/dz = 0, 
one finds the chemical potential shift 

fii = V2mlb=^. (11) 

In the last equality we have reintroduced the mesoscopic length £ and have returned to the original definition of the 
temperature gradient /3. 

The mesoscopic shift in chemical potential, = i^i: is proportional to the mesoscopic length scale £. Also note 
that fim > 0(< 0) depending on whether the (p > 0(< 0) phase is near the hot end. Physically the chemical potential 
shift (from the coexistence value /x = 0) is proportional to the temperature difference across the interfacial region, 
which is generally small. For the profile shown in Fig. 1, the chemical potential fj, evaluated from Eq. (||) is negative, 
as shown in the inset. For the </)^ model, using the parameters of Fig. 1, one finds l/z^l = 0.004/-\/2 ~ 0.00283 in 
excellent agreement with the direct simulation. At the time shown in Fig. 1, the interface hasn't quite equilibrated 
in this diffusive system. Small order parameter changes in the evolution are emphasized in the chemical potential. 
For complete equilibration, global diffusion is required. Simulations on smaller systems and in one-dimension reveal, 
more rapidly, a constant chemical potential. Recall that the homogeneous case havi ng r — tq = const has 11 — 0. 

The above analysis arriving at the leading behavior of a flat interface, i.e., (j) = ^/tJx) ta,nh{z{x)) can be viewed as 
taking the sharp interface limit in which the fast variations of the order parameter along the normal at the interface are 
considered. Note that the slow variation associated with the temperature field is effectively ignored while integrating 
though the interface. Furthermore, if the normal to the interface is tilted by an angle 6 with respect to the x— axis 
(the direction of temperature variation) one repeats the above analysis using (jj = TO(a;' cos 0) tanh(2:(a;')) where x' is 
the coordinate in the normal direction. Everything goes through as above with the mesoscopic correction becoming 
fii cos{9) where /ii is given in Eq. (pl|). Below the calculation will be slightly modified to deal with a gently curving 
interface. 

Numerical simulations, such as shown in Fig. 1, confirm that the models exhibit this shift, supporting the local 
equilibrium picture of interface in an inhomogeneous, slowly varying temperature field. As noted, a single two-phase 
interface, owing to the conservation of order parameter, can achieve an equilibrium state in a large, finite system. A 
slab consisting of a combination of two opposite interfaces, i.e., a kink-antikink, cannot. In fact the mesoscopic shift 
in chemical potential alone drives the slab to the hotter side. This phenomenon will be addressed elsewhere p9|]. 



IV. INTERFACE EQUATIONS 

We now consider a gently curved interface in the coarse-grained modeling and derive equations for the evolution 
of the interface. We restrict ourselves to two dimensions for simplicity since appropriate generalizations for higher 
dimensions can be made straightforwardly. 
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To determine the effect of an applied gradient, we first assume that locally the interface is flat and is in local 
equilibrium with x variation replaced by variation along the normal direction. The second assumption is that the 
interfacial motion arises from local quasi-stationary evolution of the order parameter 0(r, t) = (f){r — Vt). In particular, 
the evolution of the interface at any point s is characterized by the normal velocity V{s) ■ fi, where n is a unit normal 
(pointing from the 'minus' phase to the 'plus' phase). This means that order parameter or fluid flow along the 
tangential direction at any point on the interface do not contribute to the evolution of the interface. The equations 
for diffusive and fluid cases differ since the mechanisms for the interfacial motion are different. Interface equations 
in the absence of imposed gradients have been well studied in literature jl^ ; our interest here is to extend the 
derivations and applicability of such equations to inhomogeneous perturbations which can, for example, drive steady 
state behavior. 



A. Diffusive case 

Eq. (^) describing the diffusive dynamics can be inverted (see, e.g., Ref. using the assumption of local quasi- 
stationarity : 

dVG(f-r')^^^ =-Mr,<) + A(<), (12) 

where G is the Green's function for the Laplacian satisfying —'^'^G{r — r') = S{r — r'). The function A preserves the 
order parameter conservation and satisfies V^A = with zero-gradient boundary condition as does fi. One may view 
A as the change in chemical potential due to the presence of complicated interfaces in the system. The idea is to derive 
an equation for the interfacial motion starting from Eq. (|2|) by integrating the fast variation of (f> along the normal 
at the interface. A separation of length scales allows one to take sharp interface limit in which the radius of curvature 
is much greater than the interface thickness; an alternative procedure for organizing corrections has been used to 
derive interfacial dynamics in a driven diffusive system. ]20| These methods allow, however, a significant variation of 
temperature over the scale of a droplet. 

First, the local quasi-stationarity assumption allows us to replace d(j>/dt by —V ■ V0. Since the normal variation 
of (/) at the interface dominates, we approximate it further by —Vnd(j)/dn at the interface with Vn being the velocity 
along the normal. Second, we multiply Eq. (|l^) on both sides by d(f>/dn and then take the sharp interface limit. At 
any point s on the interface, by ignoring the slow tangential variation of along the interface, we can write 

^.dHdn^[-e^^-rix)^ + ^^]^-iMs)i^r , (13) 

where k is the curvature. On integration through the interface, the first term in the above is exactly the equation we 
have for the fiat interface. Hence in the sharp interface limit the first term yields the mesoscopic correction, /i,„(n-a;), 
as discussed above, while the second term gives the surface tension multiplied by the local curvature, k(s). With 
these approximations, the following integral equation results for the normal velocity at the interface 

ds'G{s\s')A{s)A{s')Vn{s') = (t(s)k(s) - /i„(n • x)A{s) - A(t)A(s), (14) 



where G{s\s') = G(r(s) — r{s')), (t{s) is the surface tension and s and s' are points parameterized by the arc length 
along the interface, /. In the above equation tangential variations in the miscibility gap A(s) give rise to to terms 
in the interface velocity which are higher order in the temperature gradient /3. The same result could be obtained 
without pausing to find explicitly the mesocsopic correction to the chemical potential. Since Eq. ( p^ will be projected 
onto the interface by action of d(t>/dn, one may use Eq. ( p^ ) with the "inner approximation" ~ 0^'' tanh(n/\/2^/), 
where the subscript / indicates the interface location. (Of course, the form of the inner approximation is justified by 
the systematic analysis leading to the mesoscopic correction.) In this paper we restrict ourselves to leading order in 
/3; hence we set A(s) — A and obtain 

ds'G{s\s')V^{s') = ^ifM^ - ("-")^: + ^(^^ (15) 
J A^ A 

Then the time-dependent Lagrange multiplier \{t) is determined as part of the solution guaranteeing conservation of 
the order parameter as specified, in the macroscopic limit, by the vanishing of the surface integral, 
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Vn{s)ds = 0. (16) 

I 

The above two equations specify the evolution of a comphcated interface in the 2D diffusive model, including the 
leading mesoscopic correction. They can be used, for example, to study the evolution of droplets as well as other 
structures and interfacial instabilities. |^l| Such equations, without the mesoscopic term, have been derived at the 
macroscopic level and used in the literature to study a variety of problems such as interfacial growth in an anisotropic 
Hele-Shaw cell ||. 

The interface equations can be used to study the evolution of complicated interfacial shapes in general. Here we 
use them to study thermocapillary response of a spherical droplet of radius R in d-dimcnsions to a small thermal 
gradient. The curvature k(s) at any point s in this case is simply given by k(s) = {d — l)/i?, and the integration 
over s' in Eq. ( p^ ) is carried over a d-dimensional sphere. Note that in the absence of a gradient the droplet remains 
stationary since, a{s) ^ <j and /i,„ = 0. As a result A is a constant and takes a value A = Aq = (j{d — l)/(i?A). This 
is immediately recognized as the required shift in the chemical potential of the system off bulk phase coexistence to 
maintain a droplet of radius R. In macroscopic treatments which assume local equilibrium, this shift is embodied 
in the the Gibbs-Thomson relation. The treatment here amounts to an alternative derivation of this relation from a 
coarse grained starting point. Other derivations are contained, for example, in Refs. 

In the presence of a thermal gradient, the surface tension a has spatial variation which drives the droplet. This 
is the standard Marangoni effect, but in a purely diffusive system. In particular, if we assume a linear variation, 
(t(s) « ctq + <7if3Rcos{9), where 9 is the angle of the outward normal to the gradient direction, and neglect higher 
order corrections in f3, there is a steady-state solution for the motion of the droplet. The velocity of the droplet 
can be obtained by solving the integral equation (p^). In particular, Vn = Vq cos(0) is a solution with A = Aq, 
corresponding to rigid motion of the droplet along the direction of the gradient. Following some of the manipulations 
in the Appendices of Ref. [Q, one expands the Green's function for two points on a sphere of radius R in spherical 
harmonics as 

oo m— -t-Z ^ 

A7r\a-a'\ + 1 

where 9, (j) locate a and similarly for a' . Inserting this in the integral equation for the drop velocity with Vn = V cos(0), 
which takes the form 

'\ d2 ,n' _ -CriPid - 1) C0S(6') ^irn C0S{9) 



~V J G{d,a')cos{9')R''dn 
one may obtain the velocity, Vq. The calculation in two dimensions is similar, and one finds 

with C = 2 in d — 2 and C = 6 in d = 3. Note that typically cti < since surface tension normally decreases with 
increasing temperature, thus drives the droplet to hotter side. It is interesting that the mesoscopic correction with fj,i 
defined in Eq. ( [ll| ) is always positive, indicating that it also drives the bubble toward the hotter side. Note that the 
dimensionality dependence only appears in the coefhcient C{d), and that the mesoscopic correction does not modify 
the scaling with droplet radius. 

As a check one can compare the results of numerical simulation of the 2D diffusive model with the above analytic 
expression. The simulations represent a direct forward integration of Eg. (pi )- (A sample plot is shown in Fig. 5a of 
Ref. [^) Good agreement with the R~^ dependence as predicted in Eq. (|lq ) is found supporting the validity of the 
interface equation approach based on local equilibrium. Including the mesoscopic correction for (3 = 0.002 one finds 
from Eq. (pj) Vb = 1.56 x IQ-" for R = 18 while the direct simulation yields Vo ~ 1.87 x 10"'', while for R = 24, 
Eq. ( p^ yields Vq = 2.8 x 10~^ to be compared with 2.71 x lO^'^ from direct simulation. It should be noted that 
the mesoscopic term in Eq. ( [l^ ) is significant. We have also been able to observe the mesoscopic correction effects in 
the numerical simulations for smaller droplets. A plot of the chemical potential itself most directly reveals this, as in 
Fig. 1. 

B. Hydrodynamic interactions 

The first assumption we make here is that the diffusive time scale of the order parameter is so large that it can 
be neglected compared to the viscous time scale set by the fluid. This is equivalent to saying that the interface is 
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being primarily advected by the fluid and that diffusive effects (necessary for full approach to equilibrium) act as 
corrections. Next, similar to the diffusive case, we appeal to quasi-stationary motion of the interface. Note that, one 
has to take both normal and tangential velocities of the fluid at the interface into account. However, the fluid velocity 
normal to the interface determines the evolution of the interface as in the diffusive case. 

The velocity equation can again be inverted, consistent with the quasi-stationary assumption, yielding the 
following equation: 



v4r) = - y I]^r,^(f,r')(AiV</))0c;V', (19) 

where Tq/j is the Oseen tensor, the form of which ensures the divergence free nature of the velocity field. This also 
implies that the normal velocity Vn on any closed surface satisfies / Vnds = 0, which, since diffusion is neglected 
as noted above, is consistent with the conservation of the order parameter content enclosed within. Note that the 
pressure term drops out in the above equation since it has been chosen to guarantee divergence free flow. In Eq. (p^), 
the /iV0 term is large mainly at the interface. One can integrate this term though the interface, similarly to the steps 
carried out for diffusive case, and arrive at the following interface equation: 

Va(s) ^^J ^0TMs')Fp{s')ds' , (20) 

where s,s' are any two points on the interface and T{s\s') = T{'r{s) — r{s')). The force -F/3(s) at the interface is 

obtained by integrating the fiVcj) through the interface. The normal and tangential forces Fn and Ft along the unit 
normals n and t at any point s on the interface are given by 

Fn{s) = -K{s)a{s) + flrriMf^ ■ (21) 

Ft(s) = ^+MmA(f.i), (22) 

OS 

where k is the curvature and a is the surface tension. Note that in higher dimensions the tangent surface is multi- 
dimensional. Since the gradient is applied along one specific direction, we can work with a single tangent vector 
by exploiting the azimuthal symmetry. The first terms on the right hand side in the above two force equations are 
macroscopic and are identical to the terms phenomenologically included by assuming a sharp interface between the two 
phases in the standard macroscopic hydrodynamic analysis. Q The terms involving are mesoscopic corrections, 
which will occur for finite size droplets. Hence the coarse-grained models reproduce the macroscopic results exactly. 
The tangential force derives macroscopically from the presence of a tangential variation of the surface tension, which 
itself follows in situations such as those with an imposed thermal gradient. 

Now we use the above interface equation (|o|) to obtain the velocity of a spherical droplet of radius i? in a small 
applied thermal gradient /3. Note that in the absence of gradient, the droplet remains stationary, and the pressure field 
satisfies Laplace's law in the model consistent with the macroscopic expectations. By assuming an approximate linear 
variation of the surface tension with temperature, a point to which we will return below, i.e., a{s) — (To+/3Rai cos{9), 
we can integrate the interface equation to obtain the droplet velocity. In the quasi-stationary approximation the center 
of mass of the droplet is assumed to move with a constant velocity V along the gradient in the lab frame, leading to 
a normal velocity V^cos(6'), where 9 is the azimuthal angle. This is due to the fact that, in the center of mass frame, 
the droplet must be stationary, and, hence, the normal velocity must vanish. Even though the component of the fiuid 
velocity is zero along the normal, the tangential velocity need not be zero in CM frame. This turns out to be the case 
for thermocapillary motion of a droplet in the hydrodynamic case. 

The CM velocity of the droplet V can be obtained by solving the integral equation in three dimensions using 
specifically the Oseen tensor Ta^(f) — g^{Saf3 + -75^). In three dimensions a useful identity |lj] is 

and for / 1, m = the integral equals 4cos(^)/15i?. Using such results one finds for the CM velocity, 

V^~-^{a,(3R) + ^{f^i^R)- (24) 

The first term is identical to the solution obtained via macroscopic analysis of the Navier-Stokes equation Q for the 
special case of two fiuids with identical fluid and thermal properties. The second term represents the mesoscopic 
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correction for the droplet velocity, where ni was evaluated in Eq. (jl^). In two dimensions the Oseen tensor is 
logarithmic, and one needs to introduce a screening length a presumably set by the droplet size. Explicitly, Tap^r) = 
— j^[5af3log{r / a) — ^^^^) . With the simplest assumption a ^ R, similar manipulations yield the center of mass velocity 
of the droplet as 



Interestingly, the macroscopic part, i.e., the first term on the right, does not depend on the screening length, a. 
For general screening length, the mesoscopic, second, term is multiplied simply by the factor (1 + ln(i?/a)). For 
a suspension of many droplets, it is reasonable to expect that a ^ R, but for a single droplet, the choice remains 
problematic. As pointed out elsewhere |^ the R dependence of the macroscopic term is expected to be dimensionality 
independent. In the coarse grained model we are considering here, the surface tension decreases with increasing 
temperature, i.e., , cti < 0. Hence the droplet is driven along the direction of the gradient towards the warmer side. 
Note that the mesoscopic term also drives the droplet to the warmer side with the scaling of velocity with radius of 
the droplet being unchanged; the total proportionality constant is shifted, however. 

As a check on our ability to analyze the motion of a droplet in the mesoscopic model, we have numerically in- 
vestigated the average velocity of the droplet and the velocity of the fluid at the interface of the droplet in the 
two-dimensional case. The simulation methodology and data arise from calculations described in Ref. |q|, in which 
Eqs. (0) and (||) were forward integrated. The macroscopic term in Eq. ( p5|) yields V ~ 0.012 for the 4^ model for 
droplet size R — 2A { and operating conditions tq = 0.5, P = 0.001, 1/ri = 5.96) which correspond to a simulation in 
that reference. The numerical agreement is excellent; Fig. 3 of Ref. Q shows V ~ 0.0005ii! ~ 0.012. This agreement 
does not address the ambiguity of the screening term in the mesoscopic correction in two dimensions. In Fig. 2, we 
show the velocity Vx and Vy at the interface of the droplet as a function of the azimuthal angle 9. One can clearly 
see the variation consistent with a non-zero tangential fluid velocity at the interface. One can show that the macro- 
scopic part of the tangential velocity is V sm(9) leading to = T^cos(20) and Vy = — Vsm{26), confirming non-solid 
body motion of the droplet in the fluid case, gj] This behavior agrees with macroscopic analysis based completely 
on the steady solution of the Navier-Stokes equation {i.e., , by redoing the analysis of Refs. for comparison 
with the two-dimensional simulations.) The order of magnitude of V obtained in the simulations is in agreement 
with the analytic estimate in Eq. ( p5| ) above. Additional quantitative investigations are left for future explorations. 
Appropriate checks should involve three-dimensional calculations to avoid the ambiguity of the screening length in 
the mesoscopic correction. The presence of the mesoscopic correction can be verified using the chemical potential 
profile which shows linearity inside the droplet along the gradient direction (i.e., V^/x is small consistent with the 
assumptions), corresponding to the mesoscopic shift in the chemical potential at the interface. 

The inclusion of a body force such as gravity into the interface equations for two fiuids having different densities 
has been discussed earlier in Ref. iQ, but without a thermal gradient. In the presence of a thermal gradient, it 
can be checked that the normal force at the interface gets appropriately modified (a term g{Ap) cos{6) with g 
the acceleration due to gravity and (Ap), the density difference enters), leading to an additional macroscopic term 
identical to the result of Young et al. [Q For example, in three dimensions the gravity coupling yields the following 
CM velocity for the droplet: 



with C = g{Ap). Note the different scaling of the macroscopic terms with droplet radius, R, indicating the interplay 
between a bulk body force and a force concentrated at the interface. With proper orientation and strength of the 
gradient relative to the gravitational force, the macroscopic terms can be tuned to cancel, leaving a small systematic 
drift toward the hot side due to the mesoscopic correction. Phase asymmetry (i.e., owing to a free energy without 
(j) ~* symmetry) also couples to a temperature gradient in a way similar to gravity modifying only the the normal 
force on the interface |q| . 



In this work we have investigated macroscopic equations for interfacial dynamics arising from two coarse-grained 
dynamical models. These equations agree with those derived from purely macroscopic considerations, which provides 
additional evidence for the utility of these coarse-grained models even away from the critical regime; furthermore the 
details of the analysis indicate clearly the assumptions required and potential sources of corrections. 




(25) 



(26) 



V. CONCLUDING REMARKS 
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More specifically, we derived interface equations starting from coarse-grained descriptions of binary systems whose 
dynamics are dominated by diffusion (Model-B) and by hydrodynamic interactions (Model-H) using a local equilibrium 
description. We allow for spatially varying potentials such as an imposed thermal gradient, as long as the spatial 
variation is slow on the scale of the interfacial width; such variation may nonetheless be macroscopically large so as 
to be significant over the scale of the relevant structures. The equations are shown to describe thermocapillary driven 
motion of a droplet quite well, yielding complete agreement with macroscopically derived results in the appropriate 
limit in which structures are large with respect to the interfacial width and the gradient is sufhcicntly weak to neglect 
all but the linear temperature coefficient of the surface tension. 

Clearly the interface equations derived are more general. The coarse-grained models naturally include non-linear 
dependence of the surface tension on the temperature. We also show how mesoscopic effects associated with the 
finite thickness of the interface do not shrink as the radius increases. The mesoscopic effects alone can play a vital 
role in driving thermocapillary migration in situations in which the curvature effects are less important (e.g., motion 
of a slab, or kink-antikink pair in a two phase system). We have added additional circumstantial evidence for the 
applicability of coarse-grained free energetics and coarse-grained hydrodynamics even away from criticality, for which 
the modeling was originally designed. The local equilibrium or "local steady state" ||2^ analysis and the interface 
equation approach are sufficiently general to allow application to a variety of other driven systems. 
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Figure Captions 



Fig.l: The order parameter profile resulting from direct numerical simulation in two dimensions is shown for (3 ~ 
0.004, To = 1 and ^ = 1 in units of thermal correlation length. One can think of the profile as a combination of 
two ramps having slopes proportional to /3. A system of length 100 is used. The inset shows the chemical potential 
evaluated from Eq. (^) at the corresponding time during the approach to equilibrium. 

Fig. 2: Simulation of droplet motion with hydrodynamic coupling, (a) The X— component of the velocity at the 
droplet interface (i.e., component along the direction of the gradient) is shown as a function of the azimuthal angle 
(6). The solid line indicates, for comparison, cos{29) behavior, (b) The y-component is shown. For comparison 
the solid line shows — sin(20) behavior. These indicate non-rigid motion for the droplet. The parameters used are 
(3 = 0.001, R = 24, To = 0.5. A 200 x 200 system was simulated in units of I. 
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